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Abstract 

About 4.4 million hadronic decays of Z bosons, recorded by the OPAL detector at LEP 
at a centre-of-mass energy of around ^/s = 91.2 GeV, are used to determine the mean 
charged particle multiplicities for the three light quark flavours. Events from primary u, 
d, and s quarks are tagged by selecting characteristic particles which carry a large fraction 
of the beam energy. The charged particle multiplicities are measured in the hemispheres 
opposite to these particles. An unfolding procedure is applied to obtain these multiplicities 
for each primary light quark flavour. This yields 

(nu) = 17.77 ± 0.51 ti'lo ' = 21-44 ± 0.63 +|-|^ , (n,) = 20.02 ± 0.13 , 

where statistical and systematic errors are given. The results for (n^) and (na) are almost 
fully statistically anti-correlated. Within the errors the result is consistent with the flavour 
independence of the strong interaction for the particle multiplicities in events from the 
light up, down, and strange quarks. 
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1 Introduction 



The flavour independence of the strong coupling is a fundamental property of quantum 
chromodynamics (QCD). A breaking of the flavour symmetry should only occur due to 
calculable mass effects. These mass effects have been observed for bottom quarks |1|,|^ 
using event shapes and jet rates in the final state of electron-positron annihilation into 
bb. Another observable which can be employed to test flavour independence is the multi- 
plicity of charged hadrons in jets originating from quarks of a specific flavour. According 
to the local parton hadron duality hypothesis (LPHD) the particle multiplicity is 
related to the gluon multiplicity inside a jet which depends on the value of the strong 
coupling constant @]. A dependence of the multiplicity of charged hadrons on the quark 
flavour has been found for heavy quarks Up to now only few measurements exist on 
the flavour independence of the strong interaction in the light quark sector (up, down, 
strange quarks) [0. Within the large uncertainties due to the limited statistics, flavour 
independence of the strong interaction is supported. 

This paper presents a new, high statistics investigation of the flavour dependence of 
the strong interaction. It is based on the mean charged multiplicity determined separately 
for events of primary up, down, and strange quarks in e"'"e~ annihilation at centre-of-mass 
energies close to the mass of the Z boson. To identify the flavour of the primary quark in 
the Z decay, the leading particle effect is exploited . This assumption of a correlation 
between the flavour of the primary quark and the type of the hadron carrying the largest 
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momentum has recently received further support by measurements of the SLD ^ and 
OPAL collaborations 0. 



The study presented in this paper uses three different selections of events with leading 
Kg and mesons, and highly energetic stable charged particles, denoted "highx^;" in 
the following. Due to the different fractions of primary up, down, and strange quarks in 
these samples a statistical decomposition of the contributions from each of the three quark 
flavours is possible. The leading Kg and K^ meson selections, even though dominantly 
stemming from primary strange quarks, yield the separation of up and down quark events 
since a leading Kg (K^) meson is rarely formed from a primary up (down) quark while 
both Kg and K^ are equally likely to be produced in a strange quark event. Together with 
the highx^; selection, containing up, down, and strange quark events in approximately 
equal proportions, properties like the mean charged particle multiplicity in events of each 
of these three light quark flavours can be determined by statistical unfolding. The mean 
charged multiplicity is determined from all long-lived charged particles in the hemisphere 
opposite to the leading particle where the two hemispheres of an event are defined by 
the plane perpendicular to the thrust axis which contains the main interaction point. 
Choosing the opposite hemisphere minimizes the bias of the measured multiplicity due to 
the high energy of the leading particle. 

Before the three selections of leading particles are presented in Section |^, the OPAL 
detector and the data samples used for this study are briefly introduced in Section ||. 
The determination of the charged particle multiplicities is discussed in Section ^. The 
investigation of various sources of systematic uncertainty is detailed in Section ^. Section |] 
presents a cross-check of the analysis using flavour fractions obtained from the data. The 
results of this analysis are presented in Section ^. 

2 The OPAL detector, data, and Monte Carlo simu- 
lation 

The OPAL detector has been described in detail elsewhere |0. The analysis presented 
here relies mainly on the reconstruction of charged particles in the large volume jet 
chamber whose performance has been presented in |]10[. A solenoidal coil surrounds 
all tracking detectors. It provides a field of 0.435 T along the beam axisQ. Tracks of 
charged particles are reconstructed with up to 159 space points. Their momenta in the 
plane transverse to the beam -RXIS, Pxy 5 CcLll be determined to a precision of cTp^y/pxy = 



y'0.022 + {p^y ■ 0.0015/GeV)^. This resolution degrades towards the acceptance boundary 
of I cos^ I ~ 0.98. From the specific energy loss dE/dx, which is measured from up to 159 
samples with a resolution of a{dE /dx)/ (dE/dx) ~ 0.035, the type of the charged particle 
can be identified over a wide momentum range. 

The analysis is based on data recorded with the OPAL detector between 1990 and 

^The coordinate system of OPAL has the z axis along the electron beam direction, the y axis points 
upwards and x towards the centre of the LEP ring. The polar angle 9 is measured with respect to the z 
axis. 
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1995 comprising about 4.4 million hadronic events at centre-of-mass energies -y/i around 
91.2 GeV (LEP I). The events considered for this study were preselected by a standard 



selection for high multiplicity events |Tl[] which relies on a minimum number of measured 
tracks of charged particles in the tracking detectors and clusters of energy deposited in 
the electromagnetic calorimeter. The remaining background of two-photon processes and 
r-pair events are estimated to be 0.07% and 0.11%, respectively. 

Tracks to be used for the reconstruction or identification of the leading particles and for 
the determination of the charged multiplicity were required to pass the quality selection 
cuts as detailed in |jl2|. These tracks had to have at least 20 hits in the jet chamber and 



a closest approach to the interaction point in the plane perpendicular to the beam axis of 
\do\ < 5 cm. The track momenta had to be between 0.1 and 65 GeV. Events with r-pairs 
are further suppressed by requiring at least seven tracks. 

The Monte Carlo simulation comprised about 8.5 million hadronic events simulated 
by version 7.4 of the JETSET program [O] which has been tuned to describe the OPAL 



data [jl4[. The generated events were passed through a detailed simulation of the OPAL 



detector [O] and processed using the same reconstruction and selection algorithm as the 



measured data. 



3 Selection of flavour enriched data samples 

Every event is divided into two hemispheres by the plane perpendicular to the thrust axis 
containing the interaction point. In each hemisphere we searched for Kg and mesons, 
and highly energetic stable charged particles, whose energy fraction xe = 2E/ ^/s is above 
a certain threshold. Due to the cuts, in particular on x^;, these three samples are enriched 
in events from up, down, and strange quarks in different proportions (see Table p. The 
Xe selection cuts were chosen to have samples of comparable size. 



3.1 Kg selection 

Kg mesons were selected via their decay into two charged pions using the procedure 
described in [|l6l. It was adapted for large momenta of the Kg mesons by dropping the 
specific do cut imposed on the two pions in the default procedure. Any two oppositely 
charged particles were combined. The invariant mass m^r^ was calculated assuming the 
pion mass for the particles. Kg candidates were required to have Im^o — m-jml < 60MeV 
using the mass of a Kg meson, m^o, given in [T^. Photon conversions were rejected by 



demanding mee > 100 MeV if the electron mass is assigned to both tracks. The impact of 
detector effects close to the acceptance boundary is reduced by restricting the polar angle 
of Kg candidates to | cos 6' | < 0.9. To enrich the sample in primary strange and down 
quarks and to suppress up, charm, and bottom events, a requirement on the scaled energy 
of the Kg candidate of xe > 0.4 was applied. Candidates with xe > 1-07 were rejected. 
This takes into account the 7 % momentum resolution for particles with the beam energy. 
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tag 




K± 


high X E 




Xe ran 




0.40-1.07 


0.50-1.07 


0.70-1.07 


Data 




= 


19359±139 


18979±138 


27909±167 


MC expectation 


A rtas 




19303± 99 


18947± 99 


27845±119 


u fraction [10^^] 


Jtag 




8.3±0.1 


21.7±0.2 


31.4±0.2 


d fraction [10^^] 


ftag 
Id 




15.6±0.2 


13.9±0.2 


27.8±0.2 


s fraction [10^^] 


ytag 




57.0±0.3 


53.1±0.3 


36.4±0.2 


c fraction [10^^] 


ytag 




14.6±0.2 


9.0±0.2 


3.0±0.1 


b fraction [10^^] 


/tag 




4.5±0.1 


2.3±0.1 


1.4±0.1 



Table 1: Number of tagged hemispheres in data and in the Monte Carlo simulation, scaled 
to the same integrated luminosity, and the flavour composition of the three samples found 
from the simulation after applying the reweighting procedure described in Section The 
errors are statistical only. 



Figure |T] (a) and (d) show the energy spectrum of the Kg candidates and the expected 
flavour composition of the events tagged by the Kg candidates. Table Ogives the number of 
tagged hemispheres in the data as well as the expected number and the flavour composition 
of the sample taken from the Monte Carlo simulation plus their statistical errors. The 
sizeable charm quark contribution is due to D mesons decaying into Kg. 



3.2 K± selection 



Charged kaons were identified using the energy loss measurements in the jet chamber . 



To effectively reject the contribution of charged pions and protons, only tracks in the 
central region of the jet chamber, | cos^ | < 0.72, were considered which have at least 130 
hits for the momentum measurement and a minimum of 100 hits for the determination 
of the energy loss. Tracks azimuthally closer than 0.5° to an anode wire plane of the jet 
chamber were not accepted to avoid biases from the significantly degraded resolution of 
both momentum and energy loss measurement close to these planes. 

Each track is assigned five dE/dx weights, Wh, calculated from the probabilities of 
the measured energy loss of the track to be consistent with that expected for a particle 
of type h, where h can be p(p), K''^, vr^, /i^ or e^. The weight is positive signed if 
di?/dx (measured) >di?/dx (expected) according to a particular hypothesis, and negative 
otherwise. K^ were selected by requiring: 

\wk\ > max{\wp\,\w^\,\w^\,\we\) 
\wk\ > 0.1 
and I^ttI < 0.1 

This yielded, according to the simulation, a composition of the sample of tagged particles 
of about 77% kaons, 13 % pions, and 9.5 % protons. The remainder is due to muons and 
charged hyperons. It was not attempted to correct the K^ selection for the pion and 
proton contributions directly. 
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OPAL 




0.6 0.7 0.8 0.9 1 0.6 0.7 0.8 0.9 1 

Xp,(high Xp cand.) Xp(high cand.)^. 



Figure 1: (a)-(c) Observed energy spectra of Kg, K^, and highx^; candidates. The solid 
line shows the Monte Carlo simulation normalized to the integrated luminosity of the data 
after applying the reweighting procedure described in Section p.4| . (d)-(f) Corresponding 
flavour fractions taken from the Monte Carlo simulation as a function of the cut on the 
minimum energy fraction, XE^raim of the tag particle. The xe cuts chosen for the flavour 
enriched samples are indicated by the vertical lines. 
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Figure |l] (b) shows the energy spectrum of the candidates. To enrich the sample in 
events of primary strange and up quarks, a cut on the minimum scaled energy of Xe > 0.5 
was applied as shown in Figure |I| (e) . Badly reconstructed tracks were rejected by requiring 
Xe < 1.07. Table |I] summarizes the number of tagged hemispheres, the expected number 
and the primary flavour composition taken from the Monte Carlo simulation. 



3.3 HighxE particle selection 

Even though the fragmentation of bottom and charm quarks is harder than that of the 
light quarks, the cascade decays of the b- and c-hadrons lead to many stable particles 
which consequently have a significantly softer energy spectrum. Thus the selection of sta- 
ble charged particles carrying a high energy fraction xe, calculated using the assumption 
that the particle is a pion, depletes an event sample of primary charm and bottom quarks. 
Besides the cuts on the scaled energy of 0.7 < xe < 1.07, additional cuts concerned the 
selection of well-reconstructed tracks which were required to have at least 130 hits used 
for the reconstruction, a polar angle of | cos 6' | < 0.72 and an azimuthal angular distance 
from the closest anode wire plane of at least 0.5°. Moreover, the Xr<p pei' degree of freedom 
of the fit of the track to the hits in the plane perpendicular to the beam axis had to be 
less than 2. This tight cut rejected badly reconstructed tracks, including those containing 
kinks due to an unrecognized decay in flight. 

The energy spectrum of the highx^; candidates and the flavour composition depending 
on the cut on the minimum energy of the highx^; particles are shown in Figures ^ (c) 
and (f). The slightly rising fraction of heavy quark events at very large xe is due to the 
remaining contribution of badly reconstructed tracks in this energy region. Such tracks 
occur with the same low probability for each primary flavour. 

In Table [l| the number of tagged hemispheres found in the data are given and compared 
with the expected number from the simulation. It should be pointed out that about 10 % 
of the highx^; hemispheres are also selected with the selection cuts. No rejection was 
applied, but it was checked that eliminating the 10 % had no significant effect on the 
result of this analysis. The small statistical correlation was taken into account. 



3.4 Checks of the flavour composition 

A crucial ingredient in this analysis is the knowledge of the flavour composition of the three 
tagged event samples. The JETSET Monte Carlo simulation using the tuned parameters 
of |n] predicts for each of the tagged particles in the region of high values of x^; a larger 
differential cross-section (1/A^) ■ (dN/dxE) than measured in the data. This discrepancy 
was accounted for by applying a flavour independent reweighting to each of the three 
tagged event samples in the simulation. The reweighting factors were derived from the 
ratio of the xe spectra of the tagged particles in the data and the simulation. They were 
parametrized by a function /{xe) = c ■ (1 — xeY whose form is motivated by the Lund 
symmetric fragmentation function [jl8|. Each tag type had its own reweighting function. 
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(a) 




Kg 


K± 


high X E 




data 
MC 


16236±127 
16185± 91 


23±5 
33±4 


116±11 
111± 8 


124±11 
123± 8 




data 
MC 


34952±187 
35075±134 




144±12 
151±10 


340±18 
350±14 


high X E 


data 
MC 


49160±222 
49245±159 






258±16 
273±12 



(b) 






K± 


high X E 




data 
MC 


3123± 56 
3118± 40 


14±4 
19±3 


24± 5 
23± 3 


lOlilO 
113± 6 


K± 


data 
MC 


8099± 90 
7938± 64 




2± 1 



197±14 
190± 9 


high xe 


data 
MC 


27909±167 
27755±119 






419±20 
407±12 



Table 2: Number of doubly tagged events when: (a) both tags belong to the lower 
Xe region, (b) at least one tag belongs to a higher xe region. The upper rows are the 
measured data. The lower rows present the reweighted simulation data (MC). Also the 
total number of tags are given for data and the simulation (third column). 



but 5 was found to be consistent for all three tags. The impact of this reweighting on the 
flavour fractions is negligible and was considered as an additional systematic uncertainty 
in Section |. 

The flavour composition and its sensitivity to changes in the fragmentation function 
can be checked with measured data. To this end the number of events with a single and 
a double tag have been counted and compared in measured and simulated data following 



the proposal of [0. Here a double-tagged event has exactly one tagged particle in each 
hemisphere, whereas single-tagged events have at least one tag. These counts depend 
on the sums and products of the flavour fractions. Due to the small number of tagging 
samples it is not possible to find a unique solution in this analysis. However, the single- 
tagged and double-tagged events can be used to cross-check the fiavour fractions obtained 
from the Monte Carlo simulation using measured data only. 

Double-tagged events were studied in different xe ranges of the tagged particles. Each 
of the three samples was divided into two xe regions. For each sample the lower xe region 
was chosen so that the impact from the choice of the parameters of the fragmentation 
function is small. The remaining part corresponds to the higher xe region which is 
particularly sensitive to the hardness of the fragmentation function. In detail, the lower 
Xe ranges were chosen for the Kg and samples to be 0.4 < xe < 0.6, and for the 
highx^; sample to be 0.6 < x^; < 0.7. The higher xe regions are 0.6 < xe < 1-07 for 
both kaon selections and 0.7 < xe < 1-07 for the highx^; sample. For both the lower and 
the higher xe regions the number of double-tagged events were determined. The events 
were classified according to whether (a) both tags belong to the lower xe region, or (b) at 
least one tag belongs to a higher xe region. Table ^ lists the double tag counts for data 
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and simulation, normalized to the integrated luminosity of the data after the reweighting. 
Good agreement between the double tag counts in data and simulation is found for both 
xe regions, as can be seen from Table 0. This can be quantified using the P^i' degree of 
freedom of the double tag counts which is 3.6/8 (5.3/8) in the higher (lower) xe regions. 
Since the double tag counts depend on the flavour fractions, the good agreement for the 
absolute number of tags and for the number of double tags leads to the conclusion that 
the estimate of the flavour fractions by the JETSET generator is reliable. 

Different multiplicities for events from primary u, d, or s quarks could imply different 
fragmentation functions for these quarks. To test the impact of different fragmentation 
functions another cross-check was conducted which considers a different reweighting factor 
/{xe) = c ■ (1 — xeY for one of the three light quarks. This corresponds to the case of 
a different fragmentation function for one quark flavour with respect to the other two. 
The flavour dependence was modelled by allowing in turn for a different value of S for one 
of the three light quarks while maintaining as a constraint in the fit of the 6 parameters 
the agreement of both the number of single tag counts and the xe spectra in data and 
simulation as before. 

When this reweighting was applied to the simulation, the light flavour fractions of 
the two kaon selections changed by less than 0.08, while for the highx^; selection the 
fractions changed by up to a factor of two. The comparison of the double tag counts 
after the flavour dependent reweighting revealed, however, a worse description of the data 
by the simulation than in the case of the flavour independent reweighting. This can be 
inferred from the P^r degree of freedom for the number of double tag counts in the 
whole Xe regions considered which increased to 27.7/7, 13.7/7, or 36.5/7, respectively for 
each sample, when allowing for a different value of 6 and, hence, a different fragmentation 
function for up, down, or strange quarks. These large values of per degree of freedom 
strongly disfavour a flavour dependent reweighting of the fragmentation functions and, 
therefore, the substantial changes of the flavour fractions associated with this reweighting. 
This gives further evidence that the flavour fractions derived from the simulation can be 
considered to be an accurate estimate of the flavour fractions in the data. We, therefore, 
use for our main result the flavour fractions obtained from JETSET without application of 
the reweighting procedure. The small effect on the flavour fractions due to the reweighting 
will yield a systematic uncertainty which is discussed in Section ^ 



4 Determination of charged particle multiplicities 

4.1 Determination of flavour dependent multiplicities 

The flavour dependent multiplicities are defined as the average number of charged particles 
with life-times r > 300 ps emerging from a decay of a Z into one of the three light 
quark flavours, u, d, or s, as for the inclusive multiplicities in [^. To obtain these 
multiplicities from the charged particle tracks measured for each of the three selections, 
several corrections must be applied. In particular, the finite detector resolution and 
acceptance, and biases due to the tag selections must be accounted for. 
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The requirement of a high energy of the tagging particles in all three selections affects 
the multiplicity in its hemisphere due to a reduced phase space for the production of 
other particles. In the opposite hemisphere the effect is much smaller. This can be 
seen in Figure § where the mean charged particle multiplicities per hemisphere in the 
three selections are plotted in bins of the tagging particle's energy. To keep the bias 
of the multiplicity due to the tagging small we determine the multiplicities only in the 
hemispheres opposite to the tagged particles. 

The remaining small dependence of the multiplicity in the opposite hemisphere on the 
energy of the tagging particle is due to the impact of the high energy of this particle. The 
higher the energy required of the leading particle the more the phase space is reduced, 
which is available for particle production, and, therefore, the mean charged multiplicity 
in the opposite hemisphere. It should also be noted from Figure |^ that the Monte Carlo 
simulation satisfactorily describes the effects in the opposite hemisphere due to the tag, 
although it tends to slightly overestimate the multiplicity in particular if the energy of the 
tagging particle is very high. This is due to a slight underestimation of the rate of badly 
reconstructed tracks in the simulation which particularly affects the hemisphere of the 
tagging particle if the measured energy fraction, x^, of such a track is close to or above 
the physical limit of 1. Since the differential cross-section is tiny for leading particles with 
xe close to 1, the effect of events selected due to badly reconstructed tracks on the mean 
charged multiplicity in the opposite hemisphere is negligible when considering the full 
range of xe used for the tagging. Any remaining discrepancies between data and Monte 
Carlo simulation in the opposite hemisphere are taken as a source of systematic error, 
which is treated in Section |^. 
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Figure 2: Mean charged multiplicity in a hemisphere as a function of the energy of the 
tagged particles. Open circles show the values obtained from the hemispheres containing 
the tagged particle, full circles correspond to the opposite hemispheres. The histogram is 
the expectation from the Monte Carlo simulation. 
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Figure 3: Charged multiplicity distributions of the three selections, measured in the 
hemispheres opposite to the tag particles. Overlaid are the expectations from the Monte 
Carlo simulation. 



Figure ^ shows the multiplicity distributions and Table ^ lists the averages of the 
distributions for the hemisphere opposite to the tag. In the Kg sample good agreement 
between data and the Monte Carlo is found, but in the events selected by a highly energetic 

or highxs tag the simulation predicts a slightly larger multiplicity. 



4.2 Unfolding procedure 

The three tag samples are composed of events with all five primary quark flavours acces- 
sible at LEP I. A statistical unfolding procedure was applied to obtain the mean charged 
particle multiplicities for each light quark flavour. The procedure simultaneously corrects 
for effects due to the selection cuts and for detector effects. 





tag 




K± 


high X E 


Data 

MC udscb 


("-Mc) — 


10.373 ±0.031 
10.323 ±0.022 


10.229 ±0.031 
10.362 ±0.021 


9.883 ±0.024 
10.031 ±0.015 


MC u 
MC d 
MC s 
MC c 
MC b 


{''^u,Mc) — 
(^d, Mc) = 
("-s.Mc) ~ 
("-c.Mc) — 


10.171 ±0.074 
10.187 ±0.056 
10.119 ±0.029 
10.853 ±0.057 
11.932 ±0.101 


10.377 ±0.045 
10.370 ±0.059 
10.242 ± 0.029 
10.684 ±0.070 
11.717 ±0.147 


10.089 ±0.027 
9.995 ±0.029 
9.921 ±0.025 
10.314 ±0.084 
11.660 ±0.132 



Table 3: Mean charged multiplicities and statistical errors of the three selections, mea- 
sured in the hemispheres opposite to the tagged particles after applying the reweighting 



procedure described in Section 3.4 
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flavour q 


mean multiplicity (rig, mc had.) 


u 


20.097 ±0.005 


d 


19.987 ±0.005 


s 


19.897 ±0.005 


c 


21.387 ±0.005 


b 


23.725 ±0.005 



Table 4: Mean generated charged multiplicities (rig^MChad.) as obtained from the Monte 
Carlo generator JETSET 7.4. Only charged particles with life-times r > 300 ps were 
considered. 



The measured mean charged multiplicities, (n**^^), in the hemisphere opposite to the 
tagged particle are expressed as 



n 



(j'=u,d,s,c,b 



('^(j.MChad.), 



{riq) , tag = K^, K=^, highxi? 



Here f^^^ is the fraction of events with primary quark flavour g in a sample denoted 
by "tag". It is listed in Table |^. The term in square brackets was also taken from the 
Monte Carlo simulation. It is the ratio of the mean charged multiplicity of the opposite 
hemisphere for events with primary quark flavour q in the "tag" sample, {i\mc)j listed 
in Table and the number for both hemispheres obtained without detector simulation, 
(^g, MChad.), which is givcu in Table ||. This ratio corrects (i) for the biases due to the 
tagging cuts, (ii) for the restriction to the hemisphere opposite to the tag particle when de- 
termining the charged multiplicity, and (iii) for effects from the selection of multihadronic 
events, the flnite resolution and acceptance of the detector. On average the correction 
factor in square brackets is 1/(2 x 1.007) for all flavours varying at most by ±3.7%, and 
where the factor of 2 is due to considering a single hemisphere only for the determination 
of the multiplicities denoted by "tag" . 

The flavour dependent mean charged multiplicities, (uq), can be determined by solving 
the equation system (p. The values of (ric) and (nb) were taken from our previous 



determination [21 



(ric) = 21.55 ±0.37 ±0.64 (2) 
(nb) = 23.16 ±0.02 ±0.45 , (3) 

where the errors are statistical and systematic, respectively. The solution of the equation 
system (|l|) using the measured multiplicities in the hemisphere opposite to the tagging 
particles from Table ^ yielded 



(riu) = 17.77 ±0.51 

{rid) = 21.44 ±0.63 

(ns) = 20.02 ±0.13 

where the errors are statistical only. Combining these, and considering the statistical 
correlation matrix: 
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u 


d 


s 


u 


1.00 


-0.90 


+0.06 


d 




1.00 


-0.42 


s 






1.00 



and using the branching ratios of Z — > uu, dd, ss from [22| one obtains (riuds) = 19.90 ± 
0.09 (stat.). This is in agreement within systematic errors with the results of |PT|,p3[ which 
were obtained using a different method. The mean charged multiphcities (n^) and (n^) 
are almost fully anti-correlated, since the statistical separation power for u and d quarks 
is due to the small but different fractions of these quarks in the Kg and the samples. 



4.3 Test of the unfolding procedure 



It was checked that the unfolding procedure described in Section ^]2| is able to determine 
the flavour dependent multiplicities. For this test the Monte Carlo data sample was 
split into two disjoint parts. The smaller, comprising about the same statistics as the 
measured data, was used in place of the measured data which were to be corrected using 
the larger part of the Monte Carlo events. The multiplicities (riu) = 20.05 ± 0.63 (stat.), 
(nd) = 19.98 ±0.81 (stat.), and (ng) = 19.89 ±0.16 (stat.) are in excellent agreement with 
the mean multiplicities obtained from the Monte Carlo generator (cp. Table §). 



5 Systematic study of the particle multiplicities 



In addition to the consistency check based on the double tag rates in Section and the 



test of the unfolding procedure in Section many sources of systematic uncertainties 



were investigated which may be subdivided into five groups, namely (i) fluctuations due 
to limited Monte Carlo statistics, (ii) the experimental precision of the multiplicity in 
charm and bottom events, (iii) impacts of detector effects and event selection and the 
dependence of the mean charged multiplicity on the choice of the hadronization model, 
(iv) variations of tagging cuts, and finally, (v) impacts on the flavour fractions, /*^^, due 
to variations of hadronization parameters of the Monte Carlo generator, the matching of 
the Monte Carlo xe spectra and hemisphere correlations to the data, and changes of the 
hadronization model. Apart from the items (i) and (ii), whose error contributions were 
estimated by error propagation, the general procedure to derive the uncertainty associated 
with the source was to repeat the complete analysis with one of the cuts or parameters 
varied. Any deviation from the result found for the standard set of cuts and parameters 
was interpreted as a systematic uncertainty. 

The respective error contributions of the uncertainties derived for the groups (i) to (v) 
were quadratically added to estimate the total systematic uncertainties. The individual 
error contributions are listed in Table ^. In the following sections the sources (iii) to (v) 
of systematic uncertainties will be discussed in more detail. 
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source of uncertainty 




A(nd) 


AK) 


data statistics 


±0.51 


±0.63 


±0.13 


MC statistics 


±0.37 


±0.48 


±0.09 


c, b multiplicity ||2T 


i 


±0.27 


±0.13 


±0.19 


detector Q 


±0.17 


±0.20 


±0.19 




XE-cut variation 
Xe < 0.9 


-0.16 
-0.06 


-0.79 
-rU.oo 

-0.06 


+0.22 
n 1 n 

+0.02 




Kg mass window 


-0.24 


+0.22 


+0.04 




Icos^Kil < 0.6, 0.8 








±0.14 


+0.20 


±0.04 


-H 


Nc3 > 80, NdE/d. > 80 








±0.09 


-0.14 


+0.02 


NdE/d. > 110 










-0.15 


+0.22 


-0.04 




\wk±\ > 0.5, Wp > 0.4, 




-0.06 




±0.16 


-0.23 


+0.04 


H 


1 cos6'highxBl < 0.6, 0.8 
A^cj > 120, 140 








±0.01 
-0.01 

< 0.01 


+0.11 
-0.16 
+0.03 
-0.01 


+0.05 
-0.03 

+0.01 














±0.01 


+0.11 


-0.03 


selection cuts total 


±0.54 
-0.36 


+0.53 
-0.88 


+0.24 
-0.12 


Alla = 250 ± 6 MeV 








±0.04 


±0.06 


±0.01 


Qo = 


= 1.9 ±0.5 GeV 










±0.03 


±0.05 


±0.03 


Oq = 0.40 ±0.03 GeV 

b = 0.52i°i° 

Ec = 0.031 ±0.011 








±0.04 

-0.06 
-1-0 1 9 

±0.06 


±0.05 

+0.22 
— n 9(1 

±0.06 


±0.01 

+0.04 
— n DQ 

±0.05 


eb = 0.0038 ±0.0010 
7s = 0.31 ±0.10 
(us/ud)/(s/d) = 0.45 ± 0.04 








±0.06 
+0.30 

— U.UO 

±0.03 


±0.11 
-0.34 

±0.06 


±0.02 

+0.06 
n n? 

— U.U ( 

±0.02 


Vd,n = 0.60 ±0.10 

= 0.40 ±0.05 
no tensor mesons 










±0.04 
-0.36 

+ U.Zu 

-0.58 


±0.07 

+0.39 
— u.oi 

+0.85 


±0.02 

-0.04 

_i_n nQ 
-|-U.Uo 

-0.14 


B{c 


^Kl + X): ±13.9%© 


B(h - 


^Kg±X): 


±10.0% 


±0.14 


+0.13 


+0.02 


B{c 


^K±±X) : ±11.5%© 


B{h - 


^K±±X) 


: ±21.6% 


+0.07 


±0.10 


±0.02 


reweighted xe spectra 








-0.35 


+0.17 


-0.03 


hemisphere correlations 








-0.01 


+0.02 


< 0.01 


hadronization total 


+0.46 
-1.05 


+ 1.25 
-0.56 


+0.11 
-0.20 


total systematic errors 


+U.86 
-1.20 


+1.46 
-1.17 


+0.39 
-0.37 


total errors 


+1.01 
-1.31 


+1.59 
-1.33 


+0.41 
-0.40 



Table 5: Compilation of individual contributions to the total error. 
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5.1 Uncertainties from detector simulation and event selection 



The detector simulation was needed to correct for the fraction of multiphcity that was not 
recorded by the detector either owing to hmited acceptance, biases of the event selection, 
or interactions in the detector such as 5-electrons or hadronic interactions, all of which 
have been studied in great detail in |^ . In total an error of about 1 % due to this 
source of systematic uncertainty was estimated in ||20|. This was also adopted for this 
measurement. 



5.2 Uncertainties from tagging 

The measured multiplicity depends on the energy of the tagged particle. Although this 
effect is reduced by measuring the multiplicity in the hemisphere opposite to the tag 
residual effects have to be taken into account. 

• Variations of the xe cuts in steps of 0.05 by as much as ±0.1 around the default 
cut value were done. The largest up- and downward excursions yielded the error 
estimates for the multiplicities. 

• Also a more stringent cut on the maximum accepted momentum, xe < 0.9, for all 
tags did not change the measured multiplicities by more than 0.4%. 

For each of the three tag types the impact of the choice of the most important selection 
cuts was studied. 



• For the Kg the allowed mass window around the world average mass value |17| was 
altered from 120 MeV within the range of 80 to 200 MeV. 

Since the selection relies very much on the capability to precisely measure the energy 
loss in the jet chamber, the relevant cuts were investigated. 

• As the precision of the measurement of dE/dx degrades towards the very forward 
and backward regions, the |cos^k±| < 0.72 requirement was varied in small steps 
between 0.6 and 0.8, taking the r.m.s. to estimate the contribution to the systematic 
uncertainty. 

• The minimum number of both jet chamber hits and dE/dx hits was relaxed from 
130 to 80. 

• The required number of dE/dx hits was increased from 100 to 110. 

• From the various modifications to the energy loss weight cuts, namely choosing a 
tighter cut on the kaon dE/dx weight of 0.5 instead of 0.1, adding an additional 
rejection cut on the proton dE/dx weight of 0.4, or considering a better pion re- 
jection by requiring its dE/dx weight to be larger than —0.06, only the first one 
contributes a significant systematic error. 
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The highx^; selection critically depends on the ability of the jet chamber to reliably 
measure high momentum particles. 

• As for K"*^, the allowed range of | cos^highx^l was varied from less than 0.72 to either 
less than 0.6 or 0.8. 

• A variation of the minimum number of hits required for a highx^; particle in the jet 
chamber by ±10 around the standard value of 130 had only a negligible effect on 
all flavours. 

• Tightening the cut on the Xrip of the track fit from 2 to 1.5 affected the multiplicity 
by less than 0.5 %. 



5.3 Uncertainties due to impacts on flavour fractions 

The Monte Carlo event generator requires several parameters to be adjusted for a proper 
description of the measured data. Changing these hadronization parameters affects the 
flavour composition, f^^^, and, therefore, tests the sensitivity of the result of this analysis 
on the particular choice of these parameters. In the following we consider those param- 
eters which are expected to have an impact on the flavour fractions. These were varied 
about their tuned values within the intervals quoted in apart from from the b and 7^ 
parameters where larger variations were investigated, in particular the range for 7^ given 
in 1^ . Most of the parameter variations did not significantly affect the flavour dependent 
multiplicities, (rig^MChad.), or the flavour fractions, f^^^, of Eq. (0). We therefore mention 
only the significant changes of the flavour fractions due to the parameter variations. 



The amount of gluon radiation was modified by varying the Alla parameter and 
the cut-off of the parton shower, Qq. 

During the hadronization step a hadron receives extra transverse momentum whose 
size is controlled by the parameter which was changed from its tuned value. 

The fraction of energy and momentum transferred from a heavy c or b quark to 
the related hadron also has an impact on the background contributions from these 
heavy quarks. The relevant e parameters of the Peterson et al. fragmentation 



function were varied. Only the change of Sc, except for the charm quark fraction 
itself, had a noticeable effect on the strange quark fraction of about —0.01 to +0.005 
for the two kaon selections. 

To account for a remaining uncertainty from the hardness of the light quark frag- 
mentation functions despite the impact of the reweighting procedure, the parameter 
b of the Lund symmetric fragmentation function in JETSET was varied in several 
steps in the range of 0.26 to 0.92. This changed in particular the strange quark 
fraction in the range of —0.06 to +0.04 for the Kg and in the range of —0.04 to 
+0.03 for the selection. The u and d fractions varied owing to this variation of 
b by ±0.012 at most in all three selections. 
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Since the statistical separation between u and d quarks is based on the Kg and 
production, the relative amount of ss quark pairs created from the vacuum, 7^, was 
changed in steps in the range of 0.21 to 0.41 |^^. This led to changes of the light 
quark fractions in the (Kg) selection ranging from —0.05, —0.03, +0.06 (—0.02, 
-0.03, +0.04) to +0.06, +0.02, -0.05 (+0.01, +0.02, -0.04), respectively, for u, d, 
s quarks. For the highx^ selection only changes of the u and s quark fractions were 
observed at a level of ±0.01. 

For baryons with strangeness an extra suppression factor (us/ud)/(s/d) exists in 
JETSET whose variation resulted in only negligible changes of both the flavour 
composition and the flavour dependent multiplicities. 

The production of vector mesons was considered. Although the down and strange 
quark fractions of the highx^; selection changed by about ±0.02, the relative portion 
Vu,d of vector mesons from u and d quarks was of minor importance for the flavour 
dependent multiplicities. 

The fraction of vector mesons from s quarks, V^, substantially affects the yields of 
leading Kg and K^ mesons. The variation of this parameter mostly changed the 
fraction from up and down quarks in the K^ selection by about ±0.01, and at the 
same level also that from down and strange quarks of the highx^; selection. Since 
the statistical u-d quark separation strongly relies on the u and d quark fractions in 
the two kaon selections, a change of this parameter led to a large error contribution 
as can be seen from Table |^. 

Similarly, the decay products of tensor mesons have a softer xe spectrum. However, 
the yields of these are not well known. Thus, the impact was estimated by analyz- 
ing about 3.5 million fully simulated events from the JETSET generator tuned to 
describe the data, but without the production of tensor mesons. This yielded a sub- 
stantial reduction of the strange quark fraction by —0.025 in the Kg, —0.035 in the 
K^, and —0.065 in the highx^; selections. The up and down quark fractions in the 
Kg and K^ selections changed by —0.01 and —0.005, whereas the d quark fraction 
in the highx^; selection increased by about +0.03. Thus the d and s quark fractions 
in the highx^; selection are about equal if tensor mesons do not contribute. The 
difference in the results with simulated tensor mesons is assigned as a systematic 
error. 

The uncertainty of the inclusive branching ratios of charm and bottom to Kg and K^ 
was considered. From the data in |17] the relative errors on the B{c KgX, K^X) 



branching ratios have been determined to be ±11.5%, ±13.9%, and on B(h 
KgX, K^X) to be ±21.6 %, ±10.0 %, respectively. Their impact on the flavour de- 
pendent multiplicities was assessed by varying the charm and bottom contributions 
in the Kg and K='= tagged event samples according to these percentages. 

To account for the poor simulation of the xe spectra of the tagged particles in the 
tuned Monte Carlo the reweighting functions of Section /{xe) = c ■ (1 — xeY, 
were employed. Applying these specific reweighting functions to all three selections 
resulted in a small reduction of the strange quark fraction by 1 to 2 % while the 



19 



heavy quark contribution is increased accordingly. Moreover, it leads to an anti- 
correlated change of the order of 1 to 3 % of the u and d quark fractions but in 
opposite directions for the Kg and the event samples. The final multiplicities of 
u and d quarks change by up to 2 %. 

• Finally, the small discrepancies between simulation and data noticed in the correla- 
tions of the opposite hemispheres shown in Figure ^were considered. Depending on 
the xe of the tagged particle the Monte Carlo events were reweighted to match the 
data. This had only a negligible effect on both the flavour fractions and the flavour 
dependent multiplicities changing them by less than 0.1%. 



In addition to the JETSET Monte Carlo event generator also the HERWIG pro- 
gram was examined. While HERWIG yields a good description of the shape of the 
Xe spectra of the tagged particles, it does not reproduce the absolute rates. For instance, 
HERWIG predicts a much larger contribution to the higha;^; sample from d quarks and 
less from u quarks due to a larger number of protons produced in HERWIG. Since the 
HERWIG expectation of the total proton yield is not in agreement with the measure- 
ments |ll7|,|27[l, one might suspect that the flavour fractions obtained from this Monte 



Carlo generator are not reliable. Thus, no systematic error contribution from using HER- 
WIG is quoted. 



6 Cross-check with flavour fractions derived from dou- 
ble tag counts 

The flavour fractions, fg'^^, can be derived from data only using the number of single 
and double tag counts since they are related to the sum and products of flavour fractions 
T9| . The non-linearity of the equation system requires for its solution more than the 



six relations which are given by the numbers of single and double tags considered in this 
analysis even when the heavy flavour contributions are kept flxed. 

Such an investigation has been done in using the single and double tag counts of 
charged pions, kaons, protons. Kg and A for xe > 0.4, 0.5, and 0.6. This analysis yielded 
the flavour dependent production rates of these hadrons where only the SU(2) isospin 
flavour symmetry and the flavour independence of the strong interaction were assumed 
for the production of charged pions, and charged and neutral kaons. 

To compare with the results from the flavour fractions, Z***^, were calculated from 
the quoted production rates neglecting correlations for this cross-check. The calculation 
took into account that these rates had been corrected for mis-identiflcation of the tagging 
particle. The highx^; sample was approximated by adding the production rates of the 
charged particles, i.e. pions, kaons, and protons. The comparison of the flavour fractions 
at an XE-cut of 0.4 for Kg, 0.5 for K^ and 0.6 in the case of the highx^; sample revealed 
some substantial differences with respect to the JETSET prediction. The Kg fraction 
from up quarks is 0.4 times and that from down quarks 1.7 times the value obtained from 
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the Monte Carlo simulation used for this analysis. The up quark contribution to the K 
sample is 1.23 times higher, while that of down quarks is about 0.75 times that obtained 
from the simulation. The fraction of high xe particles originating from up quarks is about 
10 % less than expected from the simulation. The strange quark fractions agree within 
±0.02 for each tagged event sample. These changes in the up and down quark flavour 
fractions are considerably larger than what was found from the variation of the JETSET 
parameters in Section 

Although the analysis in |p explicitly assumed flavour independence of the strong 
coupling and SU(2) isospin symmetry, the flavour fractions obtained from P] were used 
for the unfolding procedure described in Section as a cross-check of the results. This 
yielded changes of the flavour dependent multiplicities of A(?2u) = +0.60, A(nd) = —1.41, 
and A(ns) = +0.32. Since these changes are essentially covered by the systematic errors 
discussed in Section | they are not considered as an additional systematic uncertainty. 



7 Results 



Adding the individual contributions to the systematic uncertainty from Monte Carlo 
statistics, the mean charged multiplicity in c and b quark events, the correction of detec- 
tor acceptance and resolution, the selection cuts, and the parameters of the hadronization 
model in quadrature, the total systematic error of the flavour dependent mean charged 
multiplicities is about ±7 % for up and down, and ±2 % for strange quark events. The 
final results are 

= 17.77 + 0.51 (stat.)+5|[j(syst.) 
(rid) = 21.44 + 0.63 (stat.)+}-|^(syst.) 
(n^) = 20.02 + 0.13 (stat.) (syst.) . 

All errors apart from selection cuts, detector effects, and hadronization models were ob- 
tained by error propagation and, therefore, obey the statistical correlation as quoted 
in Section |4.2| . The correlation coefficients for the remaining two groups of systematic 
uncertainties were determined separately for each individual systematic variation. 

These flavour dependent mean charged multiplicities agree within the total errors, 
even though deviations may be expected due to decays of hadrons. For instance, the 
decay of the 0(1020) to charged kaons is enhanced over the decay to neutral kaons. Thus, 
the results from the JETSET generator displayed in Table || show a 1% spread between 
the charged multiplicities from up and strange quark events. 

Some of the flavour dependent mean charged multiplicities are strongly correlated 
which must be considered when comparing the results. Taking these statistical and sys- 
tematic correlations into account the ratios of the multiplicities are 

= 0.829 + 0.047 (stat.)i[j;5^J(syst.), 
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0.934 ± 0.030 (stat.)Io;o89 (sjst.) , 
1.127 ± 0.033 (stat.)lo o77 (syst.) , 




which are all close to unity. 



8 Summary and conclusions 

A measurement of the charged multiplicity in events of u, d, and s quarks from e"'"e~ 
annihilation at the Z pole is presented. Exploiting the leading particle effect to identify the 
primary quark flavour, samples of events differently enriched in u, d, and s were selected 
by tagging highly energetic Kg, K''^ and charged particles. The multiplicities per light 
quark flavour were obtained from a statistical unfolding of the multiplicities measured in 
the hemispheres opposite to the tagged particles. Hemisphere correlations were observed, 
but were found to be small and well-described by the Monte Carlo simulation. The flnal 
corrected multiplicities were determined to be 



where the (n^) and {rid) multiplicities are highly statistically anti-correlated (about —90 %). 
The ratios of pairs of these multiplicities which take the correlations into account are all 
close to unity. This agrees with the expectation of QCD that the charged particle mul- 
tiplicities are flavour independent for the light up, down, and strange quarks apart from 
small effects due to particle decays. 
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